Nodal Superconductivity with Multiple Gaps in SmFeAsOo.gFo.i 
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We report the observation of two gaps in the superconductor SmFeAsOo.gFo.i (F-SmFeAsO) with 
T c = 51.5.?f as measured by point-contact spectroscopy. Both gaps decrease with temperature and 
vanish at T c and the temperature dependence of the gaps are described by the theoretical prediction 
of the Bardeen-Cooper-Schrieffer (BCS) theory. A zero-bias conductance peak (ZBCP) was ob- 
served, indicating the presence of Andreev bound states at the surface of F-SmFeAsO. Our results 
strongly suggest an unconventional nodal superconductivity with multiple gaps in F-SmFeAsO. 

PACS numbers: 74.20.Rp, 74.45. +c, 74.50.+r 



I. INTRODUCTION 



Recently, superconductivity was discovered in Fe- 
based layered superconductor LaFeAs[Oi_ x F x ] (F- 
LaFeAsO)^. By replacing La by some other rare 
earth elements, T c was improved quickly up to 55K 2 
which is substantially higher than that of the single- 
layered cuprate superconductors and exceeds the theo- 
retical value predicted by the conventional BCS theory 3 .. 
The superconductivity was also observed in hole doped 
case^£. Thus a new family of high-T c superconduc- 
tors was opened up, providing an unique chance to un- 
derstand high-T c superconductivity. To detect the gap 
structure and pairing symmetry is an essential step to 
reveal the mechanism of these Fe-based superconduc- 
tors. The specific heat measurement on F-LaFeAsO 
showed a nonlinear magnetic field dependence of the elec- 
tronic specific heat coefficient as expected by a nodal 
superconductor—. This was proved subsequently by the 
measurements of point-contact spectroscopy^, lower crit- 
ical field H c i 8 , London penetration depth A and spin- 
lattice relaxation rate 1/Ti— The muon spin relax- 
ation measurements also presented the possibility of dirty 
d-wave pairingii. Theoretically, some calculations sup- 
port rf-wav o 12 ' 13 i 14 while others support an unconven- 
tional s-wave pairin g 15 i 16 . Most surprisingly, recent An- 
dreev reflection data favor an isotropic single gap in 
SmFeAsO0.85F0.15 and NdFeAsOo.sEr^, while a nodal 
supercoductivity with two-gap structure was suggested 
by the NMR experiment on PrFeAsO0.89F0.11—- Ac- 
tually, multi-band superconductivity has been predicted 
theoreticall y 12 ' 20 i 21 ' 22 i 23 and was supported by the mag- 
netic properties of SmFeAsOo.sFo.2 single crystal 2 ^ and 
the high-magnetic field resistance of LaFeAsO0.89F0.11—- 
However, multiple gaps have not been detected for the 
moment by both Andreev reflection or tunneling experi- 
ment which is a powerful tool to measure the supercon- 
ducting gap. Therefore, more experiments on the sam- 
ples with higher quality are strongly desired. 

In this paper, we present the point-contact spec- 
troscopy data of the compact and rigid SmFeAsOo.gFo.i 
(F-SmFeAsO) samples. A ZBCP was observed repeat- 
edly, indicating the nodal gap structure of F-SmFeAsO 
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FIG. 1: (Color online) Temperature dependence of AC sus- 
ceptibility of F-SmFeAsO. Inset: DC susceptibility of F- 
SmFeAsO. 



similar to that of F-LaFeAsO. Moreover, two different 
gaps were observed and both of them decrease with tem- 
perature and vanish at T c . These results strongly suggest 
that F-SmFeAsO has an unconventional nodal supercon- 
ductivity with multiple gaps. 



II. EXPERIMENT 

The superconducting F-SmFeAsO samples were pre- 
pared by a high pressure synthesis method 26 . The de- 
tailed information about the synthesization is elaborated 
in a recent paper—. As shown in Fig. [TJ DC susceptibil- 
ity (measured under a magnetic field of 1 Oe) and AC 
susceptibility data (measured using an AC amplitude of 
0.1 Oe) exhibit a sharp magnetic transition. The width 
defined between the 10% and 90% cuts of the transition 
is below 2K, with the middle of the Meissner transition 
at 51.5 K, indicating the good quality of the supercon- 
ducting phase. Compared with the samples synthesized 
by the common vacuum quartz tube synthesis method, 
the samples studied here are much more compact and 
rigid, thus more suitable for point-contact spectroscopy 
measurements. The point-contact junctions are prepared 
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FIG. 2: (Color online) (a) Gap structure (d-wave type) used 
to fit the experimental data, (b)-(j) Experimental data (solid 
circles) measured at various locations on the sample surface 
and theoretical calculations (solid lines) with the gap function 
shown in (a). 



by carefully driving the Pt/Ir alloy or Au tips towards 
the sample surface which is polished and cleaned before- 
hand. The tip's preparation and the details of the exper- 
imental setup were described elsewhere^?. Typical four- 
terminal and lock-in techniques were used to measure the 
conductance- voltage (dl/dV — V or G — V) characteris- 
tics. Each measurement is comprised of two successive 
cycles, to check the absence of heating-hysteresis effects. 



III. RESULTS AND DISCUSSIONS 

Figure [2] shows the G — V curves measured at various 
locations on the sample surface, which have been nor- 
malized by the high-bias data. A distinct ZBCP can be 
seen in Figs. (Hb),(e) and (f), indicating the existence 
of surface Andreev bound states, which is a clear signa- 
ture of the pair potentials with reversal sign in momen- 
tum space and is known as an unique character of nodal 
superconductors^. Such ZBCP has also been observed 
in F-LaFcAsO 7 and F-NdFeAsO^, thus nodal supercon- 
ductivity is most possibly a common property of the Fe- 
based superconductors. Blonder et al£Q have proposed a 
simplified theory for the G — V curves of an s-wave super- 



conductor/normal metal junction separated by a barrier 
of arbitrary strength. The barrier strength is parameter- 
ized by a dimcnsionlcss number Z which describes the 
crossover from metallic to ideal tunnel junction behavior 
by Z = to Z = co. Obviously, this s-wave model can 
not explain the observed ZBCP. Tanaka and Kashiwaya 31 
extended the BTK model to deal with the issue of un- 
conventional pairing symmetry. In this case, the angle 
(a) between the quasiparticle injecting direction and the 
main crystalline axis was introduced as another parame- 
ter. Moreover, the isotropic superconducting gap A was 
replaced by an anisotropic gap with d x 2_ y 2 symmetry: 
A = Aocos(29), as shown in Fig. [Ha). It was then pre- 
dicted that, for Z > the ZBCP is formed for all direc- 
tions in the a — b plane except when tunneling into the 
(100) and (010) planes. This ZBCP will be suppressed 
when the quasiparticle scattering or surface roughness 
is strong enoug h 32 i 33 . It was found that all the spectra 
shown in Fig. [2] can be described very well by this ex- 
tended BTK model. Moreover, various quasiparticle in- 
jecting angles (a) are obtained for the spectra measured 
at different positions, indicating that our measurement is 
a local detection. The broadening parameter r/Acri in 
the calculations is between 0.2 and 1, coming from some 
unclear scattering mechanisms at the interface, which is 
similar to the case of cuprate superconductors. It should 
be mentioned that the gap structure with reversal sign is 
necessary to explain our data while the details of the gap 
function can not be distinguished although the d x i_ y i- 
wave symmetry was accepted here (for example, d xy -w&ve 
symmetry 3 ^ can explain our data as well) . 

Another remarkable find in the calculations presented 
in Fig. [5] is that, the determined maximum gap Ao can 
be divided into two groups, namely, a big gap of 10.5 ± 
0.5 meV and a small gap of 3.7 ± 0.4 meV, though the 
barrier strength Z and quasiparticle injecting angle a are 
random due to the diverse configuration of the grains in 
the polycrystalline samples. The small gap can not be 
explained as degradation of the sample surface since its 
value is distributed in a narrow range far below that of 
the big gap. To our knowledge, this multi-gap feature 
is observed for the first time in transition metal-based 
high-T c superconductors. To get further insight into this 
point, we have measured the temperature dependence of 
these two gaps with distinct energy scales. 

Figure [3] shows the temperature dependence of two 
types of spectra corresponding two different gaps as men- 
tioned above. At lower temperatures, two coherence 
peaks can be seen clearly accompanied by low-energy 
depression of the quasiparticle density of states. With 
increasing temperature, the peaks are suppressed and 
smeared continuously and finally the spectra become a 
smooth featureless curve around T c . These data were 
normalized (by the data of T — 60K) in order to be com- 
pared with theoretical models, as shown in Figs. BJa) 
and (b). It was found that all the data can be fitted 
very well to the extended BTK model mentioned above. 
Fig. [3Jc) shows two spectra measured at the same loca- 
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FIG. 3: (Color online) Conductance measured for various 
temperatures from 2K to above T c . (a) large gap, (b) small 
gap. All the curves except the top one are offset downwards 
for clarity. 



tion while with different junction resistances and hence 
different barrier strengths. The similar gap value about 
10 meV was obtained from these spectra, indicating that 
the detected superconductivity does not depend on the 
junction resistance. Fig. BJd) summarizes the gap val- 
ues obtained from Figs. Ufa) and (b). The temperature 
dependence of both the large gap and the small one can 
be described by the prediction of BCS theory. The gap 
value of A =10meVleads to A /k B T c = 2.3, abit larger 
than the prediction of weak-coupling d-wave BCS theory. 
Moreover, both gaps are closed around T c reflecting the 
inter-band coupling existing in this material. 

It is interesting to note that some spectra exhibit a 
two-gap feature similar to that of MgB2. As shown in 
Fig. [He), compared with the high-temperature spectrum 
measured around T c , the low-temperature one has three 
distinct features: a ZBCP, two coherence peaks, and two 
symmetric hump at higher energy. In Fig.^f), we try to 
simulate the low-temperature spectrum after normaliza- 
tion according to the high-temperature one. In the calcu- 
lation, a simple two-component BTK model was accepted 
in which the normalized conductance (G = dl/dV) is ex- 
pressed by G — w\G\ + W2G2, where G\ and G2 are the 
conductance associated with the large gap and small gap, 
respectively, wi and W2 {wi + W2 = 1) are correspond- 
ing weights of these two gaps contributing to the total 
conductance. In order to reduce the number of fitting 
parameters, the s-wave symmetry was assumed for the 
large gap while the c? a .2_ J; 2-wave symmetry for the small 
gap. It was found that the main features in the spectrum 
can be fitted very well. Although this is a rough approx- 
imation, it captures the main physics of multiple gaps 
and nodal superconductivity. Most interestingly, the de- 
termined gap values from the calculation are llmeV and 
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FIG. 4: (Color online) (a) and (b) Temperature dependence 
of normalized spectra corresponding to Figs.[2ja) and (b), re- 
spectively. The solid lines are theoretical calculations accord- 
ing to the extended BTK model, (c) Two spectra measured 
at same position while with different junction resistance, solid 
lines are theoretical calculations, (d) Temperature dependen- 
cies of the gap values determined from fits as shown in (a) 
and (b), solid lines are the prediction of BCS theory, (e) The 
spectra measured at low temperature and around T c on same 
position of the sample surface, (f) The low-temperature spec- 
tra after normalization according to the high-temperature one 
as shown in (e), the solid line is the theoretical fit. 



4.5meV, very close to the results obtained from the spec- 
tra measured at other locations. 

Besides the sample of F-SmFeAsO which is focused on 
in this work, we have also measured the point-contact 
spectra of F-LaFeAsO^ and F-NdFeAsO. In those sam- 
ples, a more prominent ZBCP was often observed while 
some confused backgrounds can also be observed occa- 
sionally compared with the data of F-SmFeAsO, thus 
quantitative analysis is more difficult. A most possible 
explanation for this difference is that, F-SmFeAsO stud- 
ied here was prepared by high-pressure technique while 
other samples were synthesized in atmospheric pressure 
and hence more fragile and loose. Therefore, there is of- 
ten a space between adjacent grains, which can be seen 
clearly in the SEM photographs. Consequently, the tip 
is easy to penetrate through the sample surface and rests 
on a pit, leading to multiple contacts between the tip and 



4 



F-SmOFeAs Tc » 52K 
I I F-LaOFeAs Tc a 26K 
V77A F-NdOFeAs Tc • 50K 

I F-PrOFeAs Tc »45K (NMR Ref[19]) 




10 11 12 13 



FIG. 5: (Color online) Statistical chart of the Ao-values de- 
termined from the spectra measured at various positions on 
the sample surfaces. The vertical axis denotes the occurring 
times for a given gap- value specified by the horizontal-axis. 



surrounded micro-crystals. On the one hand, it increases 
the opportunities to detect strong ZBCP if a nodal pair- 
ing symmetry exists. On the other hand, it will induce 
a more confused background due to the complicated in- 
terface effect. In this case, the gap values can still be 
estimated by using a convenient method prosed by Da- 
gan et al^. The method of data analysis consists simply 
in subtracting conductances measured in an applied field 
from the zero field conductance. Since the maximum of 
the density of states at the gap value should be sensitive 
to an applied field, such subtraction will leads to a dip 
around the gap value in the subtracted spectral. Using 
this method, we have estimated Ao « 3.9meV for F- 
LaFeAsO, which is consistent with the results from both 
specific heat 6 and lower critical field 8 . 

In Fig. [51 we sum up the gap values of some Fe-based 
superconductors with different T c , which were obtained 
from our point-contact spectroscopy measurements. The 
result of PrFeAsO0.89F0.11 from NMR experiment^ 9 , is 
also presented for comparison. It was noted that for 
all these Fe-based superconductors, the superconducting 
gaps can be divided into two groups with distinct en- 



ergy scales centered at 4meV and llmeV, respectively. 
Recently, the Andreev reflection data suggested a gap 
of A « 6.7meV for F-SmFeAsO with T c = 42K and 
A fss 7meV with T c — 45. 5K. However, the isotropic s- 
wave gap was accepted in these analysis. If a d-wave gap 
was assumed, the obtained maximum gap Ao is about 
8meV, which should be ascribed to the big gap obtained 
in this work. Therefore, the high T c can be achieved now 
is most possibly dominated by the big gap. However, 
the detailed gap structure is still an open issue until the 
high-quality single crystals can be obtained. 

IV. SUMMARY 

In summary, we have studied point-contact spec- 
troscopy of the junctions built up between a normal- 
metal tip and the newly discovered Fe-based layered su- 
perconductor SmFeAsOo.gFo.i. A zero-bias conductance 
peak was observed and demonstrated to be related to 
the surface Andreev bound states. Two superconducting 
gaps with different energy scales were observed which de- 
pend on the temperature in the similar way as predicted 
by BCS theory. Our data present strong evidence that 
SmFeAsOo.gFo.i is a nodal superconductor with multiple 
gaps. 

Note added: Just before submitting this manuscript we 
have noticed a scanning tunneling spectroscopy study on 
SmFeAsOo.89^ which prepared by the same group using 
the same method as that in this work. The obtained 
spectra are in good agreement with our observations in 
Fig.CJc) and (f). 
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